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Let A be subsystem of a larger system AU B, and ip be a typical state from the subspace of the 
Hilbert space 7-Lab satisfying an energy constraint. Then pA(ip) = Tr|-0)(-01 is nearly thermal. We 
discuss how this observation is related to fast thermalization of the central region (ss A) in heavy ion 
collisions, where B represents other degrees of freedom (soft modes, hard jets, collinear particles) 
outside of A. Entanglement between the modes in A and B plays a central role: the entanglement 
entropy Sa increases rapidly in the collision. In gauge-gravity duality, Sa is related to the area of 
extremal surfaces in the bulk, which can be studied using gravitational duals. 


The usual particle scattering description of thermal¬ 
ization in heavy ion collisions (HIC) considers individual 
scattering events for each degree of freedom (i.e., par- 
ton). These scattering events are implicitly assumed to 
be incoherent - largely independent of each other. In 
the standard kinetic theory picture, multiple scatterings 
of each parton, one after another, are required to bring 
their distribution to the high entropy thermal configu¬ 
ration. By the usual reasoning, a timescale larger than 
1 fm is required in order for these multiple, successive 
scatterings to occur. 

However, nucleus by nucleus scattering in HIC is a case 
of coherent scattering of two objects, each with many in¬ 
ternal degrees of freedom. We can, in principle, describe 
the evolution of the wave functional ip of the entire two 
nucleus system from asymptotic separation at early times 
to some specific post-collision moment in time, such as 
r ~ 1 fm after the initial overlap of the two nuclei. The 
wave functional then describes a superposition of each 
possible set of interactions or scatterings of individual 
partons in the system. The bulk or macroscopic proper¬ 
ties of the remnant post-collision system are determined 
by an average over this complex superposition. Only a 
fraction of the total number of degrees of freedom are in 
the central region of the collision. The density matrix 
describing its properties is obtained through a trace over 
the remainder of the complex state. 

Recent results in the foundations of quantum statis¬ 
tical mechanics suggest that the system can thermalize 
through the spread of entanglement (i.e., superposition) 
much faster than expected from the usual (incoherent) ki¬ 
netic theory. That is, the system can evolve from a very 
atypical state (two heavy nuclei, each with large boost) to 
a more typical one (entangled superposition state) over 
a short timescale. 

Below, we discuss fast thermalization in HIC (i.e., on 
timescales of r ss 1 fm or less) in relation to the com¬ 
plex nature of the entangled superposition state, and the 
resulting growth of entanglement entropy. To be pre¬ 
cise, what requires explanation is not full thermalization 
of the central region A, but rather the weaker condition 
of isotropization so that a hydrodynamic description be¬ 


comes approximately valid. At minimum, the stress ten¬ 
sor of the matter in A must assume a diagonal form in 
its rest frame. For an overview, see [l]. 

What is new in this paper: we relate fast thermal¬ 
ization in HIC to recent results in the foundations of 
quantum statistical mechanics (the properties of typical 
states). We conjecture that entanglement between dif¬ 
ferent branches of the HIC wavefunction is the primary 
driver of fast thermalization - it is an intrinsically quan¬ 
tum effect, not a semi-classical one. This conjecture is 
supported by gauge-gravity duality: the Ryu-Takayanagi 
formula relates the area of extremal surfaces to entangle¬ 
ment entropy in HIC, implying a large increase in the 
latter during the collision. 

Quantum thermalization and properties of typical 
pure states 

Typical pure states in quantum mechanics are states 
which dominate the Hilbert measure. The ergodic the¬ 
orem proved by von Neumann H 0 implies that un¬ 
der Schrodinger evolution most systems spend almost all 
their time in typical states. Typical states are maxi¬ 
mally entangled, and the approach to equilibrium can be 
thought of in terms of the spread of entanglement. 

Consider a large system subject to a linear constraint 
R (for example, that it be in a superposition of energy 
eigenstates with the energy eigenvalues all being near 
some E*), which reduces its Hilbert space from X to a 
subspace 'Hr- Divide the system into a subsystem A , to 
be measured, and the remaining degrees of freedom which 
constitute an environment B (Fig. [I]), so 7~L = T~La 
and 

p A = Pa(iP) = Tr B \ip)(ip\, (1) 

is the density matrix which governs measurements on A 
for a given pure state ip of the whole system. The entan¬ 
glement entropy Sa is S(pa ) = — Tr hi . 

It can be shown HE] (see also EH3), using the concen¬ 
tration of measure on hyperspheres [8] (Levy’s Lemma), 
that for almost all ip € Hr, 

p A {ip) ~ Tr B (p*) = n A , 


( 2 ) 
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FIG. 1: The entire system is in a pure state ip subject to a 
constraint on total energy. Tracing over the shaded region B 
yields a density matrix pA- For typical ip (which dominate 
the set of possible states), pA is nearly thermal. 

where p* = 1 r/ dR is the equiprobable maximally mixed 
state on the restricted Hilbert space Hr ( 1 /?, is the iden¬ 
tity projection on Hr and d,R the dimensionality of Hr). 
flA = Ti's (p*) is the corresponding canonical state of 
the subsystem A. The result holds as long as cIr gU, 
where cIa and ds are the dimensionalities of the Ha and 
Hr Hilbert spaces. (Recall that these dimensionalities 
grow exponentially with the number of degrees of free¬ 
dom. The Hilbert space of an n qubit system is 2 n - 
dimensional.) In the case of an energy constraint R , 
describes a perfectly thermalized subsystem with tem¬ 
perature determined by the total energy of the system 
(i.e., a micro canonical ensemble). 

To state the theorem in [¥' more precisely, the 
(measurement-theoretic) notion of the trace-norm, is re¬ 
quired, which can be used to characterize the distance 
between two mixed states pa and Q 4 : 

II Pa ~ ^a||i = (p A - Cl A ) 2 . (3) 


entire system) apart from the maximally mixed state p* is 
exponentially small in dR, the dimension of the space Hr 
of allowed states. Conversely, for almost all pure states 
ip any small subsystem A will be found to be extremely 
close to perfectly thermalized (assuming the constraint 
R on the whole system was an energy constraint). 

As mentioned, the overwhelming dominance of typi¬ 
cal states ip is due to the geometry of high-dimensional 
Hilbert space and the resulting concentration of measure. 
It is a consequence of kinematics only - no assumptions 
have been made about the dynamics. Almost any dy¬ 
namics - i.e., choice of Hamiltonian and resulting uni¬ 
tary evolution of ip - leads to the system spending nearly 
all of its time in typical states for which the density ma¬ 
trix describing any small subsystem A is nearly thermal 
EISIS]. Typical states ip are maximally entangled (i.e., 
Sa is nearly maximal), and the approach to equilibrium 
can be thought of in terms of the spread of entanglement, 
as opposed to the more familiar non-equilibrium kinetic 
equations, which describe incoherent scattering. 

Explicit demonstrations of fast quantum thermaliza- 
tion have been obtained for broad classes initial states, 
on timescales of order the inverse temperature (i.e., av¬ 
erage energy per mode) jJUj [TT] . Indeed the conceptual 
challenge is to understand the conditions which lead to 
the more familiar slow (semi-classical) thermalization. 

Since generic pure states tend to evolve into typical 
states, any mixture of pure states is likely to evolve into 
a mixture of typical states. Hence, our analysis does 
not require any specific assumptions about whether the 
system is in a pure or mixed state. If it is in a mixture, we 
simply have (classical) probabilities of finding the system 
in one of two or more typical pure states. For simplicity, 
we can assume the system as a whole is in a pure state. 


This quantifies how easily the two states can be distin¬ 
guished by measurements, according to the identity 


IIPA-fUlli =sup|| O ||< 1 Tr( / 9 j4 0-fI A 0) , (4) 

where the supremum runs over all observables O with op¬ 
erator norm ||0|| < 1. The trace on the right-hand side of 
Q is the difference of the observable averages (O) evalu¬ 
ated on the two states pa and Qa, and therefore specifies 
the experimental accuracy necessary to distinguish these 
states in measurements of O. 

The theorem then states that (for e > 0) 


Prob 


\\PA WO - Tr s (p*) ||i > e + d A d R 1/2 


< 2 exp(—e 2 d_R/ 187 r 3 ) . 


( 5 ) 


In words: let ip be chosen randomly (according to the 
Haar measure on the Hilbert space) out of the space of 
allowed states Hr ; the probability that a measurement 
on the subsystem A only , with measurement accuracy 
given in ([5]), will be able to tell the pure state ip (of the 


Application to Heavy Ion Collisions 

In HIC we let A represent mostly soft particles in the 
central region, and let B represent all other modes (Fig. 
2). In highly central collisions at RHIC, the nucleons 
comprising each heavy nucleus pass through each other, 
losing energy due to interactions. Some of this energy 
is deposited in the central region. The rapidity of the 
original nucleons drops from about 6 to 5 [12]. Therefore, 
at most about 40 percent of the total energy in the two 
nucleus system ends up in the central region, with not all 
of it registered in detectors. In most collisions, which are 
peripheral (lower centrality), much less than 40 percent 
of the total energy is deposited in the central region. So, 
the dimensionalities dA,B of the relevant Hilbert spaces 
describing A and B obey the inequality dA -C dg. 

Thus, the conditions are appropriate for the applica¬ 
tion of the results discussed above, IF the post collision 
state ip at time r « 1 fm after the initial overlap of the 
two nuclei is sufficiently typical among states of the same 
energy. Obviously the state is not fully typical - most 
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of the energy is still in highly boosted nucleons largely 
collinear with the beamline. However, the complex pure 
state ip describing all possible combinatorial scattering 
trajectories is clearly highly entangled, and tracing over 
most of its degrees of freedom generates a large amount of 
entropy, pa is therefore a high entropy state even if it is 
not fully tliermalized. It seems plausible that pa might 
be described in macroscopic terms using some kind of 
hydrodynamic approximation [SHE] , even at times as 
early as r ss 1 fm. 

Indeed, the quantum thermalization contemplated 
here generates much more entropy in pa at early times 
than could arise from the usual picture of incoherent scat¬ 
terings of individual partons. In the usual picture, one 
must explain how degrees of freedom with typical energy 
200 to 300 MeV (given by the initial conditions used in 
hydrodynamical simulations [SHIS]) are able to assume 
a nearly thermal or at least highly entropic distribution. 
The timescale r ss 1 fm allows at most a single scatter¬ 
ing with momentum transfer less than 200 MeV for each 
particle, so there is not enough time for fine adjustment 
of energies or momenta of soft modes. 

On the other hand, superposition - quantum paral¬ 
lelism - can adjust the probability distribution over soft 
energies and momenta in the density matrix pa ■ That is, 
since all possible scatterings occur in the superposition 
state, the final probability distribution takes into account 
much finer grained effects. 

To better understand the superposition state, consider 
a parton i in one of the heavy nuclei, with momentum 
Pi. This parton will pass through the opposing nucleus, 
and (in a given scattering amplitude) might (or might 
not) scatter against any of a large number of opposing 
partons. The first scattering of i leads to a range of 
possibilities for its new momentum p^\ and after each 
subsequent scattering there is another new momentum 
p-”' ) . So a sequence of values 

{p[ 1 \p? ) ,---p[ n) }, ( 6 ) 

is associated with each parton’s trajectory. Inelastic scat¬ 
tering creates many additional particles of lower energy, 
and each of these new particles has its own trajectory. 
All possible scattering histories of all partons appear in 
the superposition state ip. When most of the degrees of 
freedom in ip are traced over, the resulting pa has very 
high entropy. 

Clearly, the amount of entanglement entropy gener¬ 
ated increases with the total energy of the collision. At 
energies below some threshold, pa presumably does not 
become approximately isotropic or hydrodynamic. Simi¬ 
larly, the strength of interactions plays an important role. 
In the limit of weak coupling we would not have signif¬ 
icant particle creation, nor multiple scatterings per par¬ 
ticle. Region A would not be highly populated with soft 
modes, whose specific state (e.g., energy-momentum dis- 



FIG. 2: The collision at « 1 fm after crossing of the two 
Lorentz-contracted nuclei. Hard and collinear modes lie in 
part B of the Hilbert space. The soft central modes lie in 
part A. Tracing over all or most of the degrees of freedom in 
B leads to a high entropy density matrix describing A. 

tribution) is correlated to that of region B. In the earlier 
discussion of “typical” states the coupling does not play 
a role because the results assume that the state ip of the 
system is typical - how and whether it becomes typical 
depends on the dynamics or strength of interactions. 

For earlier work (with a somewhat different perspec¬ 
tive) on entanglement entropy and thermalization in 
strong interactions, see fTSHTFj . 

Conclusion 

Fast thermalization in HIC might be a real-world ex¬ 
ample of a novel and intrinsically quantum mechani¬ 
cal phenomenon. Rapid growth of entanglement in the 
complex superposition of all possible scattering ampli¬ 
tudes leads to a large entropy S(pa) = — TA pa^-Pa 
even if thermalization is incomplete. This observation 
might explain the applicability of hydrodynamical mod¬ 
els (isotropization) even at early times such as r sa 1 
fm. 

It would be interesting to explore these ideas further 
in AdS models, where fast thermalization and colliding 
gravitational shocks have been studied jT9] ED]. The 
Ryu-Takayanagi (RT) formula relates the entanglement 
entropy of a gauge state A to the area of the bulk ex¬ 
tremal surface 7 a which terminates on the boundary of 
A EUE2]: 

Sa = a^u) m 

In bulk collisions which correspond to the HIC (e.g., col¬ 
lisions of gravitational waves), a horizon (black hole) is 
typically formed, which causes an increase in the area of 
the extremal surface. 
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Thus, the gauge-gravity duality implies that the entan¬ 
glement entropy of the region A increases significantly 
in HIC. That is, the thermalization (or isotropization; 
i.e., approach to equilibrium) is related to entanglement. 
See PTH25] for more discussion, including examples of 
condensed matter systems which are explicit realizations 
of the quantum thermalization phenomenon. In |251I2trj 
the distinction between von Neumann entropy (which re¬ 
mains zero in the time evolution of a pure state) and 
entanglement entropy Sa (which is sensitive to horizon 
formation in the bulk, and changes in time) is discussed 
in detail in the context of solvable boundary models such 
as free fermions. 

One can consider our paper to be simply an examina¬ 
tion of how this entanglement thermalization looks from 
the gauge/QCD side in a heavy ion collision. 

Finally, we noted that fast thermalization on 
timescales of order the inverse temperature seems to be 
generic for broad classes of quantum states and Hamilto¬ 
nians mm- If, as is sometimes conjectured, all QFTs 
have gravitational duals, this result might be understood 
in terms of horizon formation: the timescale for gravita¬ 
tional collapse is of order the light crossing time of the 
resulting black hole, which is roughly the inverse temper¬ 
ature. 
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